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Abstract. Ordinary differential equations describing the self-similar collapse of a rotating, magnetized, self- 
gravitating and isothermal filament are derived. Explicit homologous solutions are studied with special emphasis 
on the bifurcation that occurs at the magnetosonic critical point. It is shown that there is a critical value for the 
toroidal magnetic field slope at the origin above which no bifurcation occurs, the solution remains homologous, 
and below which the density and the poloidal magnetic field tend to zero at large radius (envelope) whereas the 
toroidal magnetic field and azimuthal velocity relax towards a constant value. A series of spatial profiles of density, 
velocity and magnetic field, potentially useful for comparison with numerical or observational studies, is obtained 
numerically and discussed. 
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1. Introduction 

Understanding the physics of dense cores is a crucial issue 
for star formation. In this respect the shape of the cores 
is of great interest since it depends on the physics that 
drives their evolution. 

Based on a statistical shape analysis, various studies 
have shown that the observational axis ratio distributions 
derived from a catalog of dense cores (e.g. Jijina et al. 
1999) are not compatible with a distribution of oblate 
cores and compatible or marginally compatible with a 
distribution of prolate cores (David & Verschueren 1987, 
Myers 1991, Ryden 1996, Curry 2002). The same conclu- 
sion is reached by Hartmann (2002) who shows that the 
major axis of the cores in Taurus is preferentially aligned 
with the filaments in which they are embedded. 

Interestingly Jones et al. (2001) and Goodwin et al. 
(2002) have shown recently that the apparent axis ratio 
distribution of the dense cores is compatible with the cores 
being triaxial but being more oblate than prolate. However 
the physical origin of such structures is presently unclear 
and the prolate assumption appears thus to be the sim- 
plest one which is compatible with the current observa- 
tions. 
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It is the purpose of the present study to investigate 
further the collapse of such elongated structures. 

1.1. Prolate dense cores 

Static prolate cores have been found in various theoretical 
studies. 

Curry (2000) has numerically obtained solutions of the 
Lane-Emden equations having a prolate shape. These so- 
lutions are periodic along the z-axis and reminiscent of 
the quasi periodicity of the cores sometimes observed in 
filaments (e.g. Dutrey et al. 1991). 

Curry & Stahler (2001) have shown that prolate cores 
can be solutions of the static MHD equations if the struc- 
ture of the poloidal magnetic field is such that it com- 
presses the cloud. 

Finally cores permeated by an helical magnetic field 
can also be prolate (Tomisaka 1991, Habe et al. 1991, 
Fiege & Pudritz 2000c). Fiege & Pudritz (2000c) also show 
that a large sample of cores permeated by an helical mag- 
netic field have a shape compatible with the observational 
axis ratio distribution and that the polarization hole effect 
(Fiege & Pudritz 2000d), observed towards the densest 
part of the clouds, can be nicely interpreted as a purely 
geometrical effect due to cancellation of the contribution 
of the poloidal flux in the core and the toroidal flux in the 
envelope. 
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On the other hand, prolate cores can be the result of an 
instability in magnetized filaments as first investigated by 
Chandrasekar & Fermi (1953) and recently by Nakamura 
et al. (1993) who show that a filamentary molecular clouds 
permeated by an helical field of various intensity can be in- 
stable through a combination of the gravitational and the 
magnetized sausage instability leading to the formation of 
prolate cores. More recently, Fiege & Pudritz (2000a, b) 
extend this approach to molecular clouds confined by an 
external pressure. 

1.2. Aim of the paper 

In this paper we derive self-similar solutions that can 
describe the collapse of a magnetized, rotating, isother- 
mal self-gravitating filament. We obtain four ordinary 
equations in X, the self-similar variable. These equations 
are very similar to the equations obtained by Larson 
and Penston for the collapse of a spherical non rotating 
and non magnetized cloud (Larson 1969, Penston 1969) 
and extensively studied by Hunter (1977), Shu (1977), 
Bouquet et al. (1985) and Whitworth & Summers (1985). 
Our equations admit as exact solutions, some of the ex- 
plicit solutions which are obtained by Aburihan et al. 
(2001) and by the author (Hennebelle 2001). These so- 
lutions present some important restrictions since they are 
unbounded and the magnetic poloidal pressure vanishes. 
The present study demonstrates that some of these solu- 
tions are indeed bounded, since a bifurcation arise at the 
magnetosonic point. Semi-analytical solutions that have 
a non-vanishing magnetic poloidal pressure are also ob- 
tained. 

Our goal is to understand better the gravitational con- 
traction of a filament since this process could be relevant 
in the context of prestellar prolate core collapse, and to 
obtain the spatial profiles of the fields resulting from the 
dynamical collapse of a magnetized filament. These pro- 
files could then be compared with observational data and 
guide a full numerical simulation, although the filamen- 
tary geometry used in the present study is a significant 
departure from a prolate configuration. 

The second section of the paper presents the formalism 
and discusses the physical meaning of the equations. In the 
third section, we discuss explicit solutions of these equa- 
tions and we carefully address the question of the magne- 
tosonic critical points. In the fourth section, we solve nu- 
merically the ordinary differential equations and discuss 
the physical meaning of the solutions. The fifth section 
concludes the paper. 



2. Self-similar solutions 

2.1. Reduction to ordinary differential equations 

We consider the perfect MHD equations of a self- 
gravitating gas in cylindrical coordinates. We thus ignore 
the ambipolar diffusion process. 



We make the usual self-similar reduction and set, 
the usual notations: 



X = 



P(r,t) 
$(r,t) 

Ve(r,t) 
V z (r,t) 

B e {r,t) 

B z {r,t) 



C s t' 

1 R(X) 
A^G t 2 ' 

cU{x), 

C s v^(X), 

C s v e {X), 

z 

°V 



(1) 



be(X) 



where a is a real number. The radial magnetic component 
has to vanish in order to ensure zero magnetic divergence. 

The problem considered in the present study is thus ax- 
isymmetric and independent of z but for the z-component 
of the velocity field. Also our analysis leaves aside the in- 
stability (e.g. gravitational, kink or sausage instability) 
that could develop and lead to a more complex evolution. 

If i < 0, these fields describe a contraction whereas 
they describe a rarefaction if t > 0. 

In the following, we will consider only the case t < 0, 
i.e. the core contraction until the singularity formation, 
since during the subsequent evolution, it is expected that 
a disk will form, as it is the case in the numerical study of 
Nakamura et al. (1999), that cannot be described by the 
field stated by Eq. (0). 

R, (f>,v^,vg,bg and b z are dimensionless variables and 
are expected to be of order unity. The sound speed, C s , is 
about 0.2 km s _1 for a dense core having a temperature 
of 10 K. Although in principle the self-similar solutions 
have no characteristic scales, realistic initial and boundary 
conditions limit the domain of validity. The characteristic 
spatial scale is given by r < 0.01 — 0.1 pc, the characteristic 
time scale by t < 10 3-5 years and the density is > 10 5-6 
cm~ 3 . The typical magnetic intensity is about ~ 10/iG for 

Then for higher densities, the 



a density of 
relation : 



B 
Bo 



(2) 



where 1/3 < k < 1/2 is theoretically (Mouschovias 
1976, Scott & Black 1980) and observationally inferred 
(Crutcher 1999). 

B z is the poloidal magnetic field. It tends usually to 
support the cloud against the gravitational collapse. Bg is 
the toroidal magnetic component. It has not been taken 
into account in most of the collapse studies of a magne- 
tized dense core. Such fields can be produced through the 
stretching of the poloidal magnetic field. The most im- 
portant effect of the toroidal component is that it usually 
compresses the gas radially. It can also be a support in 
the axial direction. 
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The MHD equations lead with the definitions stated 
by Eqs. (|l|) to the following equations. 
Continuity: 



-2R- XR' + ^-d x (X Vtu R) + aR = Q. 

Radial momentum conservation: 
- Xv'„ + v w v' w - V l = 

R ,/ bg „ , , b z b' z 

-r + ^- xR dx{xbe) -^r- 

Azimuthal momentum conservation: 

V^Vg 



(3) 



(4) 



XVn + V^Vg + 



X 



0. 



Axial momentum conservation: 
- a + a 2 = 0. 

Poisson equation: 

±d x (X<t>') = -R. 

Azimuthal induction equation: 
-bg - Xb' e + d x (vvjbg) + abg = 0. 

Axial induction equation: 



b z -Xb' z + ^d x {Xv^b z ) 



0. 



(6) 



(7) 



(8) 



(9) 



Eq. (Q) admits the two solutions: a = and a = 1, 
which correspond respectively to a filament static along 
the z-axis and to a filament collapsing homologously along 
this axis. The first configuration has been investigated by 
Inutsuka & Miyama (1992) in the hydrodynamical case. 
In the magnetized case, a development around X ~ 
(see Eqs. |6]) implies that b z , vg and bg cannot be all to- 
gether different from zero which is an important restric- 
tion. Moreover, it is unlikely that a cloud will collapse in 
the radial direction only and not along the z-axis. Thus in 
the following, we will consider only the second case a = 1. 

With a = 1, it is easily seen with Eqs. (||) and (Q) that, 

<// = R(v m — X). (10) 

It is also easy to see with Eqs. (|J) and (JoJ) , that the density, 
R, and the axial magnetic fields, b z , are proportional: b z = 
T Z R, where T z is a real number. 

With these two last equations, the system of Eqs. (||)- 
(^), reduces after some algebra to: 

Vrr, - X 



(uw - Xy - 1 - b 2 g /R - TIR 



V K^ + R^-X ) + % R 



X X 



X 



3_\ 

XRT 



R' = 



-vLR 



R 



X X 



-v m v 



b' a 



x(v„-xy 

-<b e _ 

Vrn - X ' 



(11) 

(12) 
(13) 
(14) 



2.2. Asymptotic behavior 

In this section we derive the asymptotic behavior of 
Eqs. (|l|)-(jl^) in the limit X — > -oo. The asymptotic 
behavior near the origin X — * will be presented and 
discussed in some details in Sect. [|. 

From Eqs. (|ll])-(|l4|), it is seen that in the limit X — > 
— oo, we have: 



R(X) ~ Ai/X, 
v m (X) ~ (A l + (b?) 2 /A 1 )ln(\X\), 



(15) 



(5) v 9 (X) ~ Ve °°exp (A, + (bfy/Ai) / \n(u)/u 2 du , 



1^1 



bg(X) ~ b^ G M-{Ai + {bT) 2 /Ai)/X). 

where A\ is a negative real number. 

Consequently at large distance, the azimuthal veloc- 
ity and the toroidal magnetic component tend to constant 
values, whereas density (and poloidal magnetic field) de- 
creases as I /X. The radial velocity tends slowly to infinity. 
This is a consequence of the infinitely long filamentary 
geometry. Let us recall that the velocity of the Larson- 
Penston solutions tends to a finite value at infinity. In 
case of strong toroidal fields, v$ and bg converge quickly 
to their asymptotic values. 

In the limit X — > — oo, we have: 
b z /bg ~ T z Aibff / X — > and the dominant forces are the 
gravitational force and the toroidal magnetic one. 



2.3. Physical discussions 
2.3.1. Physical interpretations 



The system of Eqs. (|ll|)-(|l4|) has a clear physical mean- 
ing. It is very similar to the equations derived by Larson 
and Penston (Larson 1969, Penston 1969). As for the 
Larson-Penston equations, thermal pressure (second term 
of the right-hand side of Eq. [ll]) and gravity (third term) 
are included. The solutions also include rotation (first 
term), magnetic poloidal forces (fourth term) and mag- 
netic toroidal force (fifth term). 

The system of Eqs. (|TT])- (|T^) describes the contraction 
(or expansion) of a filament with an homologous velocity 
along the z-axis. The gravitational force is zero along the 
z-axis which is consistent with an infinitely long filament. 

In a more realistic situation, the cloud geometry and 
thus the the potential should evolve with time as for the 
solutions obtained by Lin et al. (1965) for a cold and un- 
magnetized cloud. 

The velocity along the z-axis does not depend on w, 
and consequently the solutions will have no clear physical 
meaning in the limit w ~ * oo. Asymptotic behaviors in 
this limit, will then not be considered in the following. 
Because of the homologous velocity field along the z-axis, 
the validity of the solutions is also restricted to a limited 
domain not too far from the z = plan, e.g. \z\ < few \X C \ 
(X c is the position of the critical point, see next sections). 
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The magnetic poloidal field is proportional to density 
which simply means that ratio between the mass per unit 
length (m = 2ir J prdr) and the poloidal magnetic flux 
(ip = j B z rdr) is constant into the cloud. It is worth not- 
ing that m/tp — > oo when t — ► 0. This is due to the fact 
that part of the gas flows along the field line. 



It is straightforward to show that Eqs. 
quire that: 

v'o = 1/2, 



l|)-(14) re- 



whereas Eq. (nil) leads to 



(17) 



2.3.2. Magnetosonic singularities 

Eqs. ©-(H) become singular if: - X) 2 - 1 - bj/R - 
T 2 R — > 0. This is the critical point already obtained in the 
Larson-Penston equations. The first term is the square of 
the fluid velocity relatively to the similarity profile, the 
second one is the sound speed, the two last terms repre- 
sent the Alfven velocity, C a . The physical meaning is that 
when the fluid velocity relatively to the similarity profile 
reaches the fastest velocity at which waves can travel, a 
singularity occurs. It is only the solutions that pass this 
singular point that have a physical meaning. The differ- 
ence from the equations obtained by Larson and Penston 
is that since the magnetic field is included, the fastest 
waves are the fast MHD waves and travel with a velocity 
equal to: yfC* + C 2 . 

In their study, Bouquet et al. (1985) show that for any 
value of 7, the critical point reached by the homologous 
solutions is a node rather than a saddle. In the isother- 
mal case, Whitworth & Summers (1995) (see also Bouquet 
1984) were able to carry out a general study of the criti- 
cal point (not restricted to the homologous solution) and 
show that if \X C \ < 1, the critical point is a saddle and 
the solution cannot pass it, whereas in the other case, the 
critical point is a node and the solutions are physical. 



3. Explicit homologous solutions and 
magnetosonic critical point 

3.1. Explicit solutions: the homologous sub-Alfvenic 
core 

In this section, we derive explicit solutions of 
Eqs. (0), U), dH), (0). These solutions have al- 
ready been obtained by the author (Hennebelle 2001) in 
an extended formulation of self-similarity (Eq. (63) with 
T z =0). For completeness and clarity we rederive them 
here. 

We look for homologous solutions and assume the fol- 
lowing spatial dependence: 



R{X) = R , 
v„{X) = v'„ jQ X, 
v e (X) = v e>0 X, 
b z (X) = b zfi = T z R , 
b {X) = b' efi X = T <t> R X. 



1/4+Ko) 2 = R a /2 + 2(b' ) 2 /R . 



(16) 



where, as in Fiege & Pudritz (2000a), T z and are re- 
spectively the poloidal and toroidal magnetic flux to mass 
ratio. 



(18) 



In Eq. (jllj) the terms independent of X cancel out 
and lead to no new constraints, whereas the terms pro- 
portional to X 2 lead to the condition stated by Eq. (|j~8|). 
Consequently, the fields stated by Eqs. ([l6]) with the con- 
straint stated by Eq. ([l7]) constitute the asymptotic form 
of the solutions of Eqs. (|U)-@ at X -> 0. The homol- 
ogous solutions require two constraints (Eqs. |l7| and |l8|) 
and there are five parameters, which leaves three inde- 
pendent degrees of freedom (e.g. the poloidal and toroidal 
magnetic flux to mass ratio and the density). 

Eq. (|l8|) has a clear physical meaning. It is simply the 
radial momentum conservation in the pressureless limit. 
The first term of the left-hand side is the square of v'^ , 
and represents the advection term, the second one is the 
centrifugal force, the first term of the right-hand side 
is the gravitational force whereas the second one is the 
toroidal magnetic pinching. At this point there is no ex- 
plicit poloidal support since the axial magnetic component 
is uniform. However, th ere is an im plicit effect due to the 
critical point (see Sect. 3^ and 3.3). 

If b' e = and v' e = 0, then R — 1/2 and v' e = 
1/2. This is very similar to the Larson-Penston solution 
that describes a spherical cloud having uniform density 



collapsing homologously, for which Rq = 2/3 and 
2/3. 

If Rq < 1/2, small values of b' e are not allowed, 
whereas if Rq > 1/2, small values of v' e are forbidden. 
In the first case, gravity is too weak and cannot produce 
the homologous collapse by itself, the condensation is in- 
duced by the toroidal pinching. In the second case, gravity 
is too strong to produce the homologous collapse and the 
cloud must be supported by the centrifugal force. 

These solutions are valid until the critical point is 
reached and describe a sub-Alfvenic (relatively to the sim- 
ilarity profile) dense core contracting homologously. 



3.2. Critical point 

A detailed analysis of the critical point can be achieved for 
the solutions stated by Eqs. (18). This has already been in- 
vestigated for the Larson-Penston solutions by Whitworth 
& Summers (1985) in the isothermal case and by Bouquet 
et al. (1985) for any values of 7. 

For the solutions stated by Eqs. dlq), we have: 



(v w — X) 2 — 1 — bg/R — T 2 Z R = 

{\ - (b^ ) 2 /Ro)x 2 -1 - t 2 z r . 



(19) 
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Thus if —^/Rq/2 < b' g < \/Ro/2, the critical point occurs 
at: 



X r = 



1 



1/4 ~(b') 2 /iV 



(20) 



whereas if (b' g ) 2 > i?o/4, there is no singular point since 
the velocity relatively to the similarity profile is always 
smaller than the Alfven speed. 

The topological nature of the critical point can be 
studied by introducing a new parameter, s, such that: 

dX 



ds 
dR 
ds 

ds 
dvg 
ds 
dbg 
ds 
where 

N = 
D 



D. 



RN 



R 



D, 



X 
-X)N, 

VgV^ 



(21) 



= -b e N, 



X) 



D, 



X 



l 

v x~ 

X? 



X) 



TlR 
X 



3_\ 
XR) 



(22) 



l-bj/R-T 2 z R. 



The linearisation of the solutions stated by Eqs. ( |18[ ) 
near the critical point, stated by Eq. (|2(]) leads to a 
relation dY/ds = [M]Y, where Y = [X, R, v m , vg, bg] 
and where the matrix [M] is given by the value of 
d(dYi/ds)/dYj at the critical point. The matrix [M] is 
easily calculated with Eqs. ( |2l| ) and the derivative of N 
and D. We have: 



dD/dY = 



dN/dY = 



v "0 r 2 



-3R r 



-X r 



-x c ,o, 

X r 



(b c S ) 2 



(23) 
26! 



XrRr X c XrR, 



where R c , v% and b@ are the values of the fields at 
the critical point. 

The study of the eigenvalues of [M ] allows to determine 
the topological nature of the critical point. The character- 
istic polynomial P, of [AI] is: 



P(A) = A 3 A 2 - X C X 



Its roots are zero (third order) and 




0. (24) 



Ai 



-X c 1 + 2V2 



1 



A 2 = ~X C 1 - 2\/2 



(^,o) 2 ' 
i?o 

Ro 



(25) 



is eigenvalue of [M] at the third order since the 
vectors (d(dYi/ds)/dY)(X c ) are linear combinations of 
(dD/dY)(X c ) and (dN/dY)(X c ) (see Eqs. Eh. 



In the neighbourhood of the critical point, one has: 



Y(s) = E i= i i5 a J V l exp(A i s) 



(26) 



where Vj are the eigenvectors and 014 are real numbers and 
where A3, A4 and A5 are the three vanishing eigenvalues. 

Consequently, since Ai and A2 are both negative (re- 
membering that (b' g ) 2 /Ro < 1/4), the critical points are 
nodes rather than saddles. The solutions are thus able to 
cross the critical point. It does not necessarily mean that 
the solutions are always physical, since after this first crit- 
ical point, it could be possible that the solutions reach a 
second critical point through which it could be unable to 
pass. Only a numerical study seems capable to give an 
answer to this question. 

From Eqs. (Elf), we have the relations: 



dvg = 



dR = 



(t& - X C )X, 
b\ 



■dX, 



R, 



-dv v 



(27) 



■it 



-rrdvzc - ^dX, 



where v% = v e (X c ), v c m = v m (X c ) and R c = R(X C ). 

Thus, when the solutions cross the critical point, vg, 
bg and R vary according to these relations. 

When the solution crosses the critical point, D be- 
comes positive. This leads to the condition: 



,„ X 2 C + RlT 2 z - (b' e o ) 2 X 2 /R 

x? + 2R 2 ri-(b>yxi/R ■ 



(28) 



According to Eqs. ( p7p and (pa), and as pointed out 
by Whitworth & Summers (1985), there is a one dimen- 
sional space of possible solutions for X < X c which are 
compatible with a given solution for X > X c . 



3.3. Expected behaviors 



The homologous solutions stated by Eqs. ( |l6| ) and (18) 
allow us to anticipate some trends that we will recover 
numerically in the next section. 

From Eq. (|20|), it is seen that, \X C \ increases with 
r z -Ro = b z and b' g Q . From Eqs. (|l6|), it is seen that the 
value of v zn (X c ), vg(X c ) and bg(X c ) increases with \X C \, 
whereas R(X C ) and b z (X c ) do not depend on \X C \. 

Consequently, the following trends are expected: 



If the density increases, \X C \, the position of the magne- 
tosonic point, decreases (remembering that b z = T Z R). 
This is simply due to the fact that the Alfven speed 
decreases and thus the velocity of the fastest waves 
decreases as well. With Eqs. (|l8|), it is seen that the 
rotation (v' g ) increases. This means that when den- 
sity increases, the gravitational force increases as well 
and must be counterbalanced by the centrifugal force. 

If the poloidal magnetic intensity increases, \X C \ in- 
creases and the value of the radial velocity and toroidal 
magnetic field at this point increase as well. 
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Consequently, the trend is that the larger the magnetic 
poloidal intensity is, the larger \v$\ and \b$\ will be 
in the super- Alfvenic part (X < Xc). This is simply 
due to the fact that the magnetic poloidal pressure 
leads to an effective sound speed equal to: f + V 2 Z R 
(see Eq. [if]) and that , Vg and Bg are proportional 
to the sound speed. 
- If the toroidal field increases, the position of the critical 
point increases and the value of the radial velocity as 
well, which simply means that the cloud is more com- 
pressed by the toroidal pinching. As we already said, 
if (b' g o) 1 ^ 2 > -Ro/2 there is no critical point since the 
velocity of the fastest waves is always greater than the 
velocity of the fluid relatively to the similarity profile. 
Rotation must increase (Eq. |l8|) in order to counter- 
balance the toroidal pinching. 

An important parameter is the value of the radial ve- 
locity at large X. According to Eq. (|l5|), the radial ve- 
locity slowly diverges when X — > — oo but, as we already 
said, the physical meaning of the solutions is unclear in 
this limit. On the other hand, Eqs. ( |l7| ) and ( po|) give 
for v VJ (X c ) the value X c /2 < — 1 (in sound speed unit). 
Consequently, it is expected that the typical radial veloc- 
ity at few X c should be as high as few times the value of 
the sound speed. 

The smallest value of |u TO (X c )| is 1 and is reached for 
the unmagnetized solution. In that case X c — —2. Let 
us recall that for the Larson-Penston solution, X c = —3 
and v r (X c ) = — 2. The present unmagnetized solutions are 
thus less dynamical, in the inner homologous part and in 
the radial direction, than the Larson-Penston solution. 

4. Numerical results 

4.1. Method 

In this section we investigate numerically the system of 
Eqs. (|TT|)-(|l4|). These equations are four ordinary differ- 
ential equations that can be easily integrated by a fourth 
order Runge-Kutta method. 

At X = 0, the system is singular and one has to de- 
velop the solutions to first order and to start the inte- 
gration at X 0. The asymptotic development towards 
X -► is given by Eqs. © and v'^ = 1/2 (Eq. 0). 
Consequently, the number of free parameters is four (Rq , 



v' e Q , T z and b' 



s Q j and there is one more free parameter 



when the solutions cross the critical point, leading to a 
total of five free parameters. It is thus difficult to explore 
all the parameter space and we will restrict the investiga- 
tion to a limited parameter range. We start at X = — 1 -5 
and use a spatial step dX = — 10 -3 . 

The integration through the critical point requires 
some care. As we already said, the solutions for X < X c 
compatible with a solution for X > X c are not unique 
and are fixed by Eqs. (|27| ) and by the conditions stated by 
Eq. (P8|) . We proceed as follow. When the critical point 
is reached (D reverse sign), we select a value of dv^ sat- 
isfying Eq. (^8[) and give to the other fields the values 



2 3 
X/X, 





0.0 
-0.1 
-0.2 

-0.3 
-0.4 








2 3 

X/X 



Fig. 1. Solutions corresponding to Rq = 1/2, v'^ = 1/2, 
b z = (full line) 0.5 (dotted line) and 1 (dashed line) 
for b' g = O.ly/RqjA. The bifurcation occurs at X = X c 
according to Eq. (|2tj|). 



stated by Eqs. (p7[). In the following we restrict our study 
to dVxa = 0. 

We will first consider homologous cores (thus following 
Eq. ^) and then relax this assumption and consider the 
more general case of non-homologous cores. 



4.2. Homologous cores and non homologous 
super-Alfvenic envelope 

The homologous solutions stated by Eqs. ( |l6| ) and (18) 
play an important role since except for the geometry they 
are equivalent to the Larson-Penston solution that has 
been found to be in reasonable agreement with numerical 
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Fig. 2. Same as Fig. @ for b' e = 0.5^i? /4. 




Fig. 3. Same as Fig. |]for b' g = 0.9^^/4. 



simulations in previous studies (Larson 1969, Hunter 1977, 
Blottiau et al. 1988, Foster & Chevalier 1993). Moreover, 
these solutions are explicit and allow us to check carefully 
the validity of the numerical integration. 

In this section we present the numerical integration 
of Eqs. with the initial density R Q = 1/2. This 

value is important since according to Eq. (|l8|), gravity is 
exactly equal to the advection term and all values of v'g 
and b'g are allowed. 

The other parameters are: v'^ = 1/2 and in order to 
explore a large range of parameters b z = 0, 0.5 and 1 cor- 
responding to T z = 0, 1 and 2 respectively and b'g equal 
to cty/Ro/4:, whith a = 0.1,0.5 and 0.9. v'g is calculated 
according to Eqs. jisl). We integrate until X = 5X C and 
display the solutions as a function of X/X c (see Eq. po|) . 
Thus, the bifurcation always occurs at X/X c = 1. 



Before presenting the numerical results, we briefly dis- 
cuss the field intensities. 

4.2.1. Field intensity 

For orientation, we estimate the typical values of the 
poloidal field in the subsonic homologous core and the 
value of the toroidal component at the critical point as a 
function of a. With Eq. (l])-(|o|), one has: 



1 — a- 



■VT+f*Ro. 



(29) 



(30) 



If a canonical density of 10 6 cm 3 in the collapsing fila- 
ment is adopted and for R a = 1/2, T z = 1, then typically 
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in the weak, intermediate and strong toroidal field cases, 
toroidal field strengths of ~15, ~90 and ~320 /iG are as- 
sociated with poloidal field of about ~ 90/iG. 

4.2.2. Weak toroidal magnetic field 



Fig. [I] displays the results for b' g = 0.1^/i?o/4. 

For T z = 0, the solution is weakly magnetized (since 
the toroidal field is weak) and except for the geometry it is 
equivalent to the Larson-Penston solution. As pointed out 
by Whitworth & Summers (1985) this solution describes 
a strong compression wave propagating inwards into the 
cloud, the bifurcation point being the head of this com- 
pression wave. 

It is seen that for X/X c < 1 the solution is homolo- 
gous as predicted by Eqs. (|l8|). At X/X c > 1 the density 
decreases, whereas \v m \, \vg\ and \bg\ increase. The density 
decreases less rapidly for large values of b z (T z ) and as dis- 
cussed in Sect. 3.3, the velocity and the toroidal magnetic 
field increase with b z (see Eq. |l^ and ^) . 

The solutions are weakly peaked with R(0)/R(5X C ) < 
2. The velocity and toroidal magnetic fields continue 
to decrease rapidly after the critical point. We have 
ih,{X D )/v„(5X D ) ~ 0.25. 

The radial velocity is roughly 8 times larger than the 
azimuthal one. When b z = 0.5 and 1 the poloidal magnetic 
intensity dominates the toroidal magnetic intensity in the 
core (by a factor ~ 5 at X c for T z = 1, and ~ 7 for 
T 2 = 2). At X ~ 5A C , the two components are almost 
equal for T z — 1 whereas for T z — 2, b z /bg ~ 1.5. 

The radial velocity at 5X C is around 4—5 C s for low and 
intermediate poloidal magnetic intensity (T z = and 1) 
and it is much higher (~ 7) for a strong poloidal magnetic 
field (r 2 = 2). 

4.2.3. Intermediate toroidal magnetic field 

The case b' e o = Q.hyfRo/A is displayed in Fig. § The 
toroidal magnetic force in the core (2b' g X/Rq) is about 
half of the gravitational force (RqX/2). Thus the toroidal 
field is important but not dominant. 

It is seen that the bifurcation at the magnetosonic 
point is stiffer. The density decreases more rapidly 
(R(Q)/R(5X C ) < 2.5) than in the previous case. This is 
due to the fact that the cloud is compressed by the toroidal 
magnetic field. 

The velocity and the magnetic field evolution is 
much flatter in the envelope than in the core, we have 
v^(X c )/v u ,(5X c ) ~ 0.3. Thus v ru (5X c ) is even slightly 
lower than for the case with weak toroidal magnetic field 
in spite of the higher value of v m (X c ). 

It is also seen that the density variation with the 
poloidal magnetic intensity is small compared to the vari- 
ation of the velocity and the toroidal magnetic component 
; we discuss this further in the next section. 

The amplitude of the radial velocity is about 2/3 the 
amplitude of the azimuthal one. The poloidal magnetic 



R - 


- R, 


V^J - 




v e - 


-> Kvg, 


b z - 


+ Kb z , 


be - 


■* Kb*. 



component dominates the toroidal one in the core (ratio 
0.5-0.75 at X = X c ). The two component become com- 
parable at X ~ 1.5X C which is earlier than for the weak 
toroidal magnetic field. At X/X c = 5 b z /be ~ 3. 

4.2.4. Strong toroidal magnetic field 

The case V 6fi = 0.9^i? /4 is displayed in Fig. § The 
toroidal force is larger than the gravitational force (factor 
~ 1.62) and the solution describes the inward propagation 
of a strong magnetic compression wave into the cloud. 

The most striking features of the solutions with strong 
toroidal magnetic field is that the density is almost in- 
dependent on the poloidal magnetic intensity, whereas 
the other fields are roughly proportional to T z . This is 
due to the fact, that since the toroidal field is strong 
(bg(X c )/R > 3), the magnetic and centrifugal forces are 
large compared to the thermal pressure. Thus in Eq. (O), 
the thermal pressure can be neglected (term — 1 in D and 
l/X in N). In this limit, it is then seen that Eq. ( [TT| ) is 
invariant under the rescaling: 



(31) 



Another interesting feature is that the radial velocity 
and the toroidal magnetic field are almost constant just 
after the critical point. They then decrease significantly 
again around X/X c ~ 1.5. The radial velocity then de- 
creases according to Eqs. (|lj]) whereas ve and be decrease 
more slowly and rapidly converge toward a constant value. 

The solutions are even more peaked than in the pre- 
vious cases, with R(0)/R(5X C ) ~ 5 whereas the radial 
velocity is flatter with v^(X c )/v zo (5X c ) ~ 0.5. 

The azimuthal velocity is about 1.3 — 1.5 larger than 
the radial velocity. The toroidal magnetic field rapidly 
dominates the poloidal one since at X/X c ~ 0.5, the two 
field intensities are comparable. At X/X c , bg/b z ~ 2 and 
at X/X c — 5, we have bg/b z ~ 10. 

For T z = or T z = 1, v^(bX c ) ~ 4 - 5 C 3 . 

4.2.5. Varying the density: forbidden values 

In this section, we vary the density of the homologous core 
in order to explore the dependence of the fields on this 
parameter. We look carefully at two particular cases and 
briefly comment on other values. In order to make compar- 
ison with previous cases easier we choose T z — 1 (dotted 
line of previous section) and be = 0.1y / i?o/4, 0.5-\/i?o/4 
and 0.9-y/ 'Rq/ '4, i.e weak, intermediate and strong mag- 
netic toroidal intensity. 

The two cases, i?o = 1/4 and Rq = 1 are more closely 
considered. In the first case, gravity alone cannot explain 
the collapse (see Eq. |l^) . The collapse is thus assisted by 
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Fig. 4. Solutions corresponding to Rq = 0.25 and b Zi o — Fig. 5. Solutions corresponding to Rq = 1 and b z $ = 1 
0.25 (r 2 = 1) for b' go = 0.5^0/4 (dotted line) and b' eo = (T z = 1) for b' g = 0.1y/R /4 (full line), b' go = 0.5^/ R / '4 
0.9a/%/4 (clashed line). (dotted line) and b' g = 0.9 v /i? /4 (dashed line). 



the toroidal magnetic field. In the second case, gravity is 
too strong and must be partially counterbalanced by the 
centrifugal force. 

The results are respectively displayed in Fig. |] and 
Fig. |. The full line is b e = O.ly^Ji (this case is not 
possible for Rq = 1/4 according to Eqs. [l8|), the dotted 
line be = 0.5^i? /4 and the dashed line bg = 0.9-y/ Rq/ 4. 

From a comparison between Fig. |4| and it is seen that 
the density is more peaked whereas the radial velocity is 
flatter for small densities (Rq — 1/4) than for large one 

(Ro = 1). 

For Ro = 1/4, the azimuthal velocity vanishes 
according to Eq. fll8[ ) when b g = O.b^Ro/A. When, 
b g = 0.9-v/ i?o/4, it is seen that the radial velocity and the 



toroidal field increase slightly just after the critical point, 
then reduce at X/X c ~ 1.5. 

Smaller values of density (Rq < 1/4), require values 
of b' g o closer to -J Ry/A in order to have positive values 
of \v$\ (see Eqs. ^). Similar behaviors to the previously 
presented cases are obtained. 

For larger values of density (Rq > 1), we find that with 
b' g = 0, the solutions oscillate around the homologous so- 
lution, i.e. after the critical point the density decreases and 
increases alternately. When b' g ^ (e.g. b' g > 0.1) then 
the solutions cross another critical point and we stop the 
numerical integration. This behavior is somehow similar to 
the results obtained by Whitworth & Summers (1985) who 
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£1 




1 2 3 4 5 

X/X, 



Fig. 6. Solutions corresponding to Rq = 0.25, 6 2i o = 0.25, 
b' o = 0.9^0/4. Full line is Sv' e = -0.05, dotted line is 
Sv' = —0.15 and dashed line Sv' = —0.25 (see text). 




-2.0 E , , , , 3 

1 2 3 4 5 

X/X 

Fig. 7. Solutions corresponding to Rq = 0.5, 6 Zj o = 0.5, 
b' g o = 0.9^i?o/4. Full line is Sv' e = -0.02, dotted line is 
Sv' g = —0.06 and dashed line 5v' g = —0.1 (see text). 



find that the allowed values of the density at the origin 
in the Larson-Penston equations are quantized. It could 
thus be possible that larger values of Rq are allowed in 
Eqs. (|Il|)-(|l^). Only a systematic numerical search could 
answer this question. 

4.3. Non-homologous cores 

In this section, we investigate the case of non homolo- 
gous cores, i.e. we start with initial conditions that do not 
satisfy the second conditions stated by Eqs. jlq). Let us 
introduce the parameter, 5v' e — v' e — v' e , where v' e is the 
value stated by Eqs. (|l^). The smaller \5v' e \ is, the closer 
we are to the homologous solutions. 



We consider the two cases Rq — 0.25 ; and Rq = 
0.5 with r z = 1 and b' e = 0.9^i? /4 for three val- 
ues of Sv' e (-0.05, -0.15, -0.25 for R = 0.25 and 
-0.02, -0.06, -0.12 for Rq = 0.5). The results are dis- 
played in Figs. H and ^. Full lines represent the small- 
est value of \Sv' g \, dotted line, the intermediate value, and 
dashed line the largest. 

It is seen that the cores are no longer homologous and 
become less flat when \6v' g \ increases. Also, the transition 
between the sub-Alfvenic part and the super- Alfvenic one, 
is smoother. The poloidal force does not vanish any more 
in the cores and acts explicitly everywhere. 

For R =0.25, the density and radial velocity values at 
5X C do not vary significantly with 5v' g , whereas the values 
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of the azimuthal velocity and the toroidal magnetic field 
decrease significantly when \8v' e \ increases. 

For larger values of \Sv' m \ the solutions (for Rq = 0.25 
and i?o = 0.5) cross another critical point. We consider 
them as not physical and we stop the numerical integra- 
tion. It is then not possible for these two cases to obtain 
solutions with a vanishing azimuthal velocity as suggested 
byEq. @. 

For larger values of Rq, the largest allowed value of 
Sv'vj decreases rapidly (for R = 1, we find that the largest 
allowed \Sv' g \ is between 0.03 and 0.04). This is consistent 
with the fact that the hom ologou s solutions with large Rq 
are not physical (see Sect. 4.2.5 ). 

It is worth noting that the value of the radial velocity 
at 5X C is around 2.5 to 3 which is lower than the values 
obtained previously. The smallest value of v^{bX c ) we 
found is around 2 and is obtained for R < 1/8, be ~ 
v/iVi, T z < 1 and \v' 9 Q \ < 0.2. 

This value is still 3 to 4 times higher than the value ob- 
served by Tafalla et al. (1999) in the starless cores L1544. 
However, it is known from previous study (e.g Foster 
& Chevalier 1993) that the Larson-Penston solution (in 
which the velocity tends to 3.3C S at large r) is in good 
agreement with the inner collapsing part of the cloud, the 
external part being significantly different. This is due to 
the fact that the self-similar solutions make the assump- 
tion that the system is infinite and is self-similar every- 
where. In particular, the self-similar assumption does not 
take into account the external medium. Self-similar solu- 
tions thus appear in the inner part of the cloud when the 
boundary conditions have been forgotten by the system. 



5. Conclusion 

In this paper, we have derived self-similar solutions able 
to describe the inner part of a collapsing, rotating, magne- 
tized, self-gravitating, isothermal filaments. The filaments 
is collapsing homologously along the z-axis, the slope of 
the axial velocity being two times the slope of the radial 
one at the origin. The four ordinary equations obtained 
are similar to the equations derived by Penston (1969) 
and Larson (1969). They present a critical magnetosonic 
point that induces a bifurcation. 

For some of the homologous solutions already obtained 
previously, (Aburihan et al. 2001, Hennebelle 2001) we 
carried out a study of the critical point and found that 
there is a critical value of b' (slope of the toroidal field at 
the origin) above which no bifurcation occurs. By studying 
the eigenvalues of the linearized equations in the neigh- 
bourhood of the critical point, we demonstrate that these 
solutions are able to cross the singularity since it is a node 
rather than a saddle point. 

We have then explored the system numerically and 
obtained a series of density, velocity and magnetic fields 
profiles that could help to understand the observational 
data and could be used as benchmarks for a full numerical 
simulation. 



The solutions obtained in this paper have some restric- 
tions. First, they are valid not too far from the equatorial 
plane since they have an homologous axial velocity field 
that quickly diverges and not too far (few \X C \) from the 
z-axis since the axial component of the velocity field does 
not depend on w. 

Nevertheless, these solutions are the first semi- 
analytical solutions describing the condensation of a mag- 
netized rotating filament and having complex (non ho- 
mologous) spatial profiles. They can be used for future 
analytical or numerical studies of the gravo-magnetic con- 
densation. 
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